Universal scaling behavior of pseudogap with doping in high Tc cuprates; temperature 

and doping dependence of spectral intensity 
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Based on our improved SU(2) slave-boson approach (Phys. Rev. B 64, 052501(2001)) to the 
t-J Hamiltonian, we report a scaling behavior of pseudogap with doping and the temperature and 
doping dependence of spectral functions. In addition we discuss the cause of hump and quasi-particle 
peak in the observed spectral functions of high Tc cuprates. It is demonstrated that the sharpening 
of the observed quasi-particle peak below Tc is attributed to the bose condensation of holon pair. 
From the computed ratios of pseudogap Ao to both the superconducting temperature Tc and the 
pseudogap temperature T* as a function of hole doping concentration x, we find that there exists a 
universal scaling of these ratios with doping, that is, the hyperbolic scaling behavior of the former, 
~ x~°' with a ~ 2 and near doping independence of the latter, -^^j- ~ 4 ~ 6 are found. 
PACS numbers : 74.25. Jb, 79.60.-i, 74.20.Mn 



Angle resolved photoemission spectroscopy(ARPES) 
measurements of high Tc cuprates revealed both the tem- 
perature and doping dependence of both spectral peak in- 
tensity and pseudogap at varying pseudogap temperature 
T* and superconducting temperature Tc. These mea- 
surements have shown a continuously increasing trend 
of pseudogap (leading edge gap) with decreasing tem- 
perature even below and appearance of sharp peaks 
below Tc [0J^. Earlier, Wen and Lee ||] reported the 
momentum dependence of the spectral function based 
on their SU(2) slave boson theory involving single-holon 
bose condensation. Chubukov et al. obtained the peak- 
dip-hump feature of the spectral function by using their 
spin fermion model [Q. Most recently, Muthukumar et 
al. studied the momentum dependence of the spectral 
function in the resonating- valence-bond state [^. How- 
ever there exists lack of comprehensive investigations on 
the temperature and doping dependence of spectral func- 
tions. Lately, based on our newly improved SU(2) slave- 
boson approach of introducing holon-pair channels |^ 
over a previous one we were able to successfully obtain 
the arch-shaped superconducting temperature line in the 
phase diagram of high Tc cuprates. Further application 
of this theory resulted in the peak-dip-hump structure 
of optical conductivity in good agreement with observa- 
tions 1^ . In the present study using this theory we report 
temperature and doping dependence of spectral functions 
for both the normal and superconducting states and dis- 
cuss the origin of the hump and the quasiparticle peak 
observed in ARPES Q Finally we report our theo- 
retical finding of the scaling behavior of the spin gap vs 
the superconducting and pseudogap temperatures with 
respect to the hole doping concentration. 

In the slave-boson representation j|] ||l^], the electron 
operator of spin a, can be written as a composite 
of spinon and holon operators; = b^/o- in the U(l) 
representation and Cq = -^h^tpa in the SU(2) theory 



with a = 1,2, where faib) is the spinon(holon) anni 

hilation operator in the U(l) theory, and 4'i — f '{f 

V J 2 
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are the doublets of 



and/.= (J^ 

spinon and holon annihilation operators respectively in 
the SU(2) theory. Judging from the presence of the in- 
tersite charge-charge coupling (jiinj) in addition to the 
spin-spin coupling (S,j • S^) in Heisenberg term of the t- 
J Hamiltonian and the consideration of on-site charge 
fluctuations it is obvious that the charge degree 

of freedom has to appear in the slave-boson represen- 
tation of the Heisenberg term. The Heisenberg term in 
the t-J Hamiltonian is, then, derived to be P(Si • Sj — 

in the U(l) theory, and -i(l - hlh,){l - h]h,){fljl^ - 

flfh){h3hr - /2,7l^) in thc SU(2) theory §. Here P 
represents the projection onto the singly occupied site 
or empty site. The four holon (boson) operator in the 
above expression represents the charge degree of freedom 
and the four fermion operator, the spin degree of free- 
dom. In the present SU(2) study the phase fluctuation 
effects of the spinon pairing order parameter are taken 
into account. 

Introducing Hubbard Stratonovich transformations for 
direct, exchange and pairing channels and a subsequent 
saddle point approximation, the t-J Hamiltonian is de- 
composed into the spinon sector, and the holon sec- 
tor, i. 
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where x-y fljfai + ,j(^il^yi {b\jbu - ^i^zj) > is the 
hopping order parameter, Af^ =< /ij/2j - /2j/ii >, the 
spinon pairing order parameter, A^^.^^ =< biabpj >, the 
holon pairing order parameter, /x^, the effective chemical 
potential, and cc, the hole doping concentration. With 
the uniform hopping order parameter, Xij — Xi the d- 
wave spinon pairing order parameter, A{^ = ±A^ with 
the sign +(— ) for the nearest neighbor link parallel to 
X (y) and the s-wave holon pairing order parameter, 
^\j.ap = '^b{Sa.iSf3.i-Sa,2S[3.2), the quasiparticlc energy 
for spinon is given by 0] 



f _ 



where the spinon single particle energy is given by, 
f_ J{\-xf 



x(cos kx + cos ky), 

and the spinon pairing gap or simply the spin gap, 
Aj = J(l — a;)^A/(cos k^ — cos ky). 



(2) 



(3) 



(4) 



The single particle (electron) propagator of interest is 
given by a convolution integral of spinon and holon prop- 
agators in the momentum space 1^ , 



dk duj 
"(2^ 



(5) 



The mean field Green's functions are G'^^^^^(k, cj) = 



< 



r[^„,(x,t)^^,„(0,0)] > and 
< r[6z(x,t)6L(0,0)] > 
respectively. The symbol < > refers to the finite tem- 
perature ensemble average of an observable quantity O, 
<0>= ^tr{e-^"0). 

The one electron removal spectral function, /(k, w) is 
obtained from 11131, 



/(k,c.) 



1 



lmG(k,uj + iO+)f{uj) 



(6) 



where /(w) is the Fermi distribution function. The 
Heisenberg coupling constant of J = 0.2 t and the hop- 
ping strength of t = 0.44 eV are chosen in the present 
calculations. Using the SU(2) theory, the predicted val- 
ues of optimal hole doping Xo, pseudogap temperature 



T* and bose condensation temperature Tc are Xo = 0.13, 
T* = 0.029i(148K) and = 0.021t (107.2if) respec- 
tively. The spectral function is obtained from the con- 
volution integral of the holon and spinon Green's func- 
tions G(k, (jj) . For the evaluation of the spectral functions 
/(k, cij), the line width of the Lorentzian function is set 
to be e = 0.01t(4.4mey) with the choice of t = 0A4eV. 

It is known from the ARPES measurements of LSCO 
[ p^ and BSCCO that a good Fermi nesting appears 
in a region near the M point (k = (tt, 0)) at which maxi- 
mal antiferromagnetic correlations between electrons are 
realized. The predicted binding energy monotonically de- 
creases with the increase of momentum k along F — M, 
in agreement with the ARPES [0. In Fig.l, predicted 
spectral functions at k = (7r/2,7r/2) below and above 
Tc are displayed. The peak did not disappear even at a 
temperature above Tc (e.g., T = 0.022t(112.3/C) with a 
choice of t = 0.44ey) and showed a decreasing trend of 
peak intensity with increasing temperature, in agreement 
with the ARPES measurement of Shen et al ||l|. At all 
temperatures below T* the zero gap is predicted. 

In Fig. 2 we display the computed spectral functions 
/(k, w) at k = (tTjO) both below and above T^. The 
computed gap size at k = (tt, 0) is predicted to coin- 
cide with the spin gap size A^ in Eq.(^. The peak 
position occurs at smaller binding energies at temper- 
atures (0.022^(112. SiC)) above T^. Although not shown 
here, a gap begins to open at higher critical temperature 
T* = 0.029t{U7.9K). This is now taken to be the spin 
gap temperature. The predicted peak position shift, that 
is, the variation of leading edge gap is much slower be- 
low Tc and reaches a finite value at T = 0, in complete 
agreement with observations ||l^ ]l7[ |. 

In order to see only the role of the antiferromagnetic 
fluctuations of the shortest possible correlation length, 
namely the spin singlet pair excitations in the evaluation 
of the convolution integral, we deliberately removed the 
contribution of the holon-pair channels which comes from 
momenta k = (0,0) and k — (7r,7r). The sharp peak 
is now seen to disappear and only the hump structure 
persists while the pseudogap remains unchanged, indi- 
cating that the pseudogap is caused by the formation of 
spin singlet pairs as is shown in the inset of Fig. 2. This 
indicates that the presence of holon-pair bosons is essen- 
tial for yielding the observed sharp quasiparticlc peaks 
at the M point below Tc- That is, the quasiparticlc peak 
is caused by bose condensation below Tc. The appear- 
ance of the hump in the absence of bose condensation is 
seen to be caused by the antiferromagnetic spin fluctua- 
tions of the shortest possible correlation length, that is, 
the spin singlet pair excitations. In Fig. 3 the doping 
dependence of spectral functions is displayed for under- 
doped, optimally doped and overdoped cases at a low 
temperature below T^ T = 0.004i (T = 20AK), near 
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the M point (k = (O.Stt, 0)) for comparison with obser- 
vation. In agreement with the ARPES |^ the pre- 
dicted spectral weight of the sharp quasiparticle peaks is 
seen to increase as the hole concentration increases upto 
a tested value in the overdoped region. The leading edge 
gap is shown to decrease with doping in agreement with 
observations At all doping rates both the peak and 
hump structures are revealed with the absence of dip. 
It is observed to be caused by electron-phonon coupling 
[ p^ , which is not considered in our present treatment of 
the t-J Hamiltonian. Figs. 4 displays as a function of dop- 
ing rate the ratios of the pseudogap 2Ao(fc = (7r,0)) at 
a temperature close to OK to the superconducting tem- 
perature Tc and the spin gap temperature T*, that is, 
2Ao vs Tc and 2Ao vs T*. It is possible that each high 
Tc cuprate(e.g., LSCO, YBCO and BSCCO) may have 
an effectively different J value causing variation in Tc- 
It is quite encouraging to find that in agreement with 
observations , ^^o(k--{7T,o)) jjgf,]-ga,ses rapidly, showing 
a hyperbolic behavior with doping concentration while 
"^^oik-i^fi)) g]-^Q^g nearly doping independence of rang- 
ing between 4 and 6, as shown in Fig. 4. It is found that 
these ratios do not appreciably change with the varia- 
tion of Heisenberg coupling J. This is due to the nature 
that Ao, T* and Tc are found to be proportional to J, 
keeping these ratios unchanged ||]. Thus we found that 
independent of Heisenberg exchange coupling there ex- 
ists a universal scaling behavior of the pseudogap with 
doping for the ratios of Aq to Tc and Aq to T*. 

Both the SU(2) and U(l) slave-boson theories pre- 
dicts sharp quasiparticle peaks at k = (tt, 0) and k = 
(7r/2, 7r/2) particularly below Tc. Although not shown 
here, the predicted leading edge gap at k = (tt, 0) showed 
a continuous increase as temperature drops down from a 
pseudogap temperature T* to temperatures below the su- 
perconducting temperature Tc. This indicates that the 
origin of the observed leading edge gap in the supercon- 
ducting state is the same as the ones in the pseudogap 
phase. That is, the pseudogap is caused by the forma- 
tion of spin singlet pairs(spinon pairs). The presence of 
the singlet pairs, namely the antiferromagnetic fluctu- 
ations of a shortest possible range is shown to persist 
even in the superconducting state. We infer from our 
computed results that the appearance of the distinctively 
sharp quasiparticle peaks below Tc is attributed to the 
enhanced probability of spin singlet pair excitations as a 
result of coupling to holon-pair bosons in the supercon- 
ducting state. It is of note that such coupling between 
the two degrees of freedom is well manifested in the last 
term of Eq.(0). 

In summary the present study revealed numerous 
salient features. They are; 1. regarding the hole dop- 
ing dependence of spectral functions near k — (tt, 0), the 
spectral intensity of quasiparticle peak shows an increas- 
ing trend with hole concentrations in the underdoped re- 



gion, 2. regarding the temperature dependence of pseu- 
dogap(spin gap) at k = (7r,0), the gap size continuously 
increases as temperature decreases from the pseudogap 
temperature T* to temperatures even below the super- 
conducting temperature Tc, 3. regarding the tempera- 
ture dependence of peak intensity, a decreasing trend of 
spectral peak intensity is observed as temperature in- 
creases, 4. the zero leading edge gap at k = (7r/2,7r/2) 
is unchanged both below and above Tc, 5. regarding the 
momentum dependence of binding energy, it decreases 
with increasing momentum k along the T — M direction, 
6. the sharp quasiparticle peak below Tc is attributed to 
the bose condensation of holon pairs, and the hump is 
caused by the antiferromagnetic spin fluctuations of the 
shortest possible correlation length, 7. the coexistence 
of the spin singlet pair excitations and superconductivity 
below Tc is predicted, and finally, 8. in complete agree- 
ment with observations, independent of the Heisenberg 
exchange coupling there exists a universal scaling behav- 
ior of the pseudogap with doping, namely the hyperbolic 
scaling behavior in and near doping independence 
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FIG. 1. Spectral functions /(k,a;) for k = (7r/2,7r/2) at the 
predicted optimal doping of Xo = 0.13 (Tc = 0.02W(107.2A'), 
T* = 0.029t{U7.9K)) with J = 0.2t at T = 0.002t(10.2/s:) be- 
low Tc and T = 0.022t(112.3is:) above Tc. 
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FIG. 2. Spectral functions I(k,Lo) for k = (7r,0) at the 
same predicted optimal doping Xo, critical temperature Tc and 
pseudogap temperature T* at the same temperatures below 
Tc and above Tc as shown in Fig.l. The inset is I{k,Lo) for 
k = (7r,0) with the exclusion of only k = (0,0) and (7r,7r) 
at T = 0.002i(10.2i4:) below Tc. 
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FIG. 3. Doping dependence of spectral functions I(k,ij) near 
the M point (k = (0.87r,0)) at temperature 0.004t(20.4iC) be- 
low Tc for underdoped(x = 0.04 and x = 0.085), optimal 
doping(xo = 0.13) and overdoped(x = 0.2) rates. 



4 



(a) 



i 


Experiment 




J=0.2t 




J=0.3t 




J=0.4t 






\ 

\. 











10 


* 


8 




6 








o 


4 


< 




M 






2 








0.5 1 1.5 
(b) 



Experiment 
J=0.2t 
J=0.3t 
J=0.4t 



0.5 1 1.5 2 

x/x 

FIG. 4. The ratios of the spinon pairing gap at T = O.OOlt 
(~ 5K) with the choice of t = 0.4461^(13] to (a) the supercon- 
ducting temperature, '^^<'^^=('^'°)) and (b) the spin gap tem- 
perature, ^^"^''^^r" as a function of where Xo = 0.13 
is the predicted optimal doping concentration for J/t — 0.2 
and 0.3, and Xo ~ 0.14 for J/t = 0.4. The experimentally ob- 
tained universal ratios for BiiSriCaCu^Os+x, La2-xSrxCu04 
and ThBttiCuOe are denoted as solid circles[20]. 
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